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PRIORITY 

This application claims priority to a U.S. Provisional Application filed on April 10, 
2000 having U.S. Provisional Application Serial No. 60/196,080; the contents of which are 
incorporated herein by reference. 
10 BACKGROUND OF THE INVENTION 

1. Technical Field 

The present invention relates to the study of neurons including neuronal 
development and the effects of various agents on neurons through analysis of optical imagery. 



15 2. Description of Related Art 

Recent large-scale genome sequencing and expression analysis have uncovered 
a multitude of genes implicated in brain development, learning and memory, regeneration, and 
neurological diseases. Determining the function of these genes and other substances 
necessitates advanced techniques for controlling the timing and location of gene expression, 

20 combined with specific assays of neuronal cell function or morphology. 

Neurons are known to include dendrites, cell bodies and axons. The areas 
between adjacent neurons are known as synapses. In most synapses presynaptic axons 
terminate on dendritic spines. Spines are small bulbous compartments consisting of a spine 
head attached via a thin neck to the dendritic shaft. Spine necks diffusionally isolate spine 

25 heads from their parent dendrites, e.g., spines restrict the diffusion of Ca 2+ and other second 
messengers. In addition, spines contain a variety of organelles. Spines contain post-synaptic 
densities, one of the most complex signaling assemblies, which include synaptic receptors and 



their regulators as well as various structural and adhesion molecules. Spines also contain 
machinery required for protein translation. 

Spine growth is associated with synaptogenesis. During periods of 
synaptogenesis dendrites grow filopodia, relatively long actin rich protrusions that often make 
5 several synapses; later these filopodia are typically replaced by mature spines. Molecular 

mechanisms underlying spine genesis and stabilization are beginning to be investigated. Since 
actin is highly enriched in dendritic spines, most studies have focused on pathways associated 
with the regulation of actin dynamics. Calcineurin appears to be important in actin stabilization 
in spines. The Rho family of small GTPases, including Rho, Rac, and Cdc42, regulates various 

10 aspects of the actin cytoskeleton and also modulates dendritic structure and spine density. 

Recent experiments have revealed that important aspects of cognitive function, 
such as experience-dependent plasticity, neural integration and learning -and-rnemory are 
correlated with variations in dendritic branching morphology, and with spine density and 
distribution. Similarly, age-related deficits in short-term memory, important forms of neural 

15 dysfunction have been localized, in part, to dendrites and spines. 

Recently direct measurements in mammalian brain slices have revealed 
that synaptic plasticity can manifest itself in sprouting of filopodia, and spines, in an N-methyl- 
D-aspartate receptor (NMDA-R) dependent manner. On the other hand prolonged NMDA-R 
activation leads to a loss of spines. Other studies have shown that spine density and shape is 

20 controlled by background electrical and synaptic activity. For example, brief exposure to the 
sodium channel blocker TTX increases the density of dendritic spines. Spines display subtle 
actin-based motility that appears to be abolished by low levels of alpha-amino-3-hydroxy-5- 
methyl-4-isoxalone (AMP A) or glutamate. Low levels of AMP A also block lesion-induced 
spine degeneration. Thus, spines are stabilized by low levels of activation of synaptic 

25 receptors, but grow in response to a global reduction of activity, as well as strong focal increase 
in activation of synaptic receptors. The signaling mechanisms underlying this complex 
response are not understood. 

These and other findings have motivated extensive efforts to obtain quantitative 
descriptions of dendritic and spine morphologies, both statically and dynamically. 

30 Due to its superior resolution capability in revealing ultrastructures at synaptic junctions, serial 
section electron microscopy (SSEM) has been used to quantify dendritic spine structures in 
three-dimensions (3-D). This is, however, a non- vital form of observation and an extremely 
labor-intensive histological approach requiring the physical sectioning of the tissue into very 
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thin sections and detailed manual and/or semi-automatic registration and outlining of the 
structures on each serial section. 

Modern fluorescence microscope methods, such as confocal laser scanning 
microscopy (CLSM) and two-photon excitation laser scanning microscopy (2PLSM) offer 
5 many advantages over SSEM, at the expense of reduced resolution. Sectioning is achieved by 
limiting the detection (CLSM) or excitation (2PLSM) of fluorescence to a sub-femtoliter focal 
volume. Optical imaging is rapid and noninvasive. The exquisite selectivity of fluorescence 
allows the detection of even single molecules against a background of billions of others. 
Optical microscopy thus occupies a unique niche in biology due to its ability to perform 
10 observations in intact, living tissue at relatively high resolution. The properties of fluorescence 
microscopy images are well understood. To image neuronal structure, neurons are labeled with 

■TiiorpcopntjronlprMilpc ihz*t fill fhp r^/tnnlacm Virimncrpnpmiclv Vrwpl valnp.Q r^nnrt thp 

convolution of the density of fluorescent probes with the point-spread function (PSF) of the 
imaging system, which is essentially equal to the focal volume and is easily measured. Studies 

1 5 of morphological plasticity based on CLSM and 2PLSM measurements of spine length and 
density have been described. 

Despite these advances in modern imaging techniques, the analysis of 
neuronal structure has remained largely manual. The considerable amount of time and effort 
required for manually extracting spine measurements has precluded routine studies of large 

20 amounts of data. In addition, results are not precisely reproducible as accuracy is dependent on 
the skill and habituation of the user. A few detection and estimation techniques (Rusakov et 
al., Quantification of dendritic spine populations using image analysis and a tilting dissector, J. 
Neurosci. Methods, 1995; 60: 1 1-21; Watzel, et al., Detection of Dendritic spines in 3- 
dimensional images, DAGM-Symposium Bielefeld, 1995; 160-167; Herzog, et al., Restoration 

25 of three-dimensional quasi-binary images from confocal microscopy and its application to 
dendritic trees, Cogswell CJ, Conchello J. and Wilson, T., editors, Three-Dimensional 
Microscopy: Image Acquisition and Processing IV. SPIE Proceedings, 1997; Kilborn et al., 
Delineating and tracking hippocampal dendritic spine plasticity using neural network analysis 
of two-photon microscopy, Soc. Neurosci. Abstr., 1988; 24: 422-425) of varying degrees of 

30 automation have been suggested to overcome the tedium and improve accuracy and 

reproducibility of the result, none of which has apparently been used and verified on large data 
sets. Rusakov, et al., supra, applied a medial axis construction (skeletonization) to 2-D 
dendritic images to obtain spine length measurements in 2-D and estimated the corresponding 

3 



3-D measurements using a stereological sampling procedure. As the medial axis is sensitive to 
surface features, manual screening of the medial axis was required to select among spine, 
dendrite and irregular surface-induced features (i.e. artifacts). Since measurements were based 
solely on the medial axis, no volumetric estimates were obtainable. Watzel, et al., supra, have 
5 also suggested the detection of dendritic spines using medial axis based identification. Their 
3-D algorithm was restricted to images containing a single dendrite. The dendrite backbone 
('centerline') was extracted from the medial axis and the remaining medial axis 'spurs' 
branching off the backbone were used to identify candidate spines. A length tolerance was 
employed to distinguish true spines from artifact 'spurs'. No further analysis beyond that for a 
10 single dendritic image was presented. Herzog, et al., supra, employed a 3-D reconstruction 
technique using a parametric model of cylinders with hemispherical ends to 'fit' the shape of 
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detect and had to be manually added to the model. An approach using neural network 
recognition for spines (Kilborn, et al., supra) has also been suggested. 
1 5 The need to understand the structure and function of neurons is a continuing 

one. Similarly, there is a continuing need to determine the effects of various substances on the 
development, structure and function of neurons. 



SUMMARY OF THE INVENTION 

20 An algorithm for determining neuronal structure by analyzing a microscopy image is 

provided wherein the algorithm includes a processing module for processing the image and 
extracting neuronal structures therefrom based on geometrical features of the neuronal 
structures, and an analyzing module for analyzing the extracted neuronal structures to 
determine at least one characteristic thereof. 

25 Also provided is a method for determining the effect of a substance on a neuron which 

includes subjecting the neuron to the substance, imaging the neuron to generate at least one 
image, subjecting the at least one image to an algorithm which contains (i) a processing module 
for processing the image and extracting neuronal structures therefrom based on geometrical 
features of the neuronal structures and (ii) an analyzing module for analyzing the extracted 

30 neuronal structures to determine at least one characteristic thereof, and comparing the at least 
one characteristic to a corresponding at least one characteristic of a control neuron. In addition, 
dual-color images can be created wherein one color is used to determine the structure of 
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dendrites and spines and another color is used to measure the distribution of various cellular 
proteins. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 FIGS. l(a)-(c) illustrate a raw microscopy image, and deblurred microscopy images 

after 5 and 20 iterations; 

FIGS. 2(a)-(b) illustrate a medial axis from a segmented image of a dendritic image 
with arrows indicating loop (L) formed by overlapping spines, separate skeletons for 
disconnected features (D) and spurious cell debris, and two backbones extracted from this 
10 medial axis; 

FIGS. 3(a)-(c) are diagrams illustrating a projected view of 12 spine candidates (shaded 
gray) with db(S)=\2 9 a sketch indicating distances and spine candidates used in the spine 
detection algorithm of the present invention; and an ideal spine candidate symmetric along the 
line SE; 

15 FIG. 4 is an illustration in 2-D of the orientation criterion used to determine if detached 

(D) or attached (A) spine components need to be merged; 

FIGS. 5(a)-(b) provide a comparison of manual and automatic spine detection on a 
segment of a hippocampal CA1 dendrite, respectively; 

FIG. 6 provides a comparison of manual and automatic measurements of individual 
20 spine length (left), average spine length (center) and spine density (right) for the dendrite 
shown by FIGS. 5(a)-(b); 

FIG. 7 illustrates charts of spine volume-length scatterplots according to determined 
spine type for all spines in experiments Ej and Ei\ 

FIGS. 8(a)-(b) illustrate charts showing the number of spines detected in the image at 
25 each time step and the detection history of 52 spines followed for 25 minutes, e.g., 1 spine was 
detected in each image taken over the 25-minute period, where spine 14 was seen sporadically 
over the entire period; 

FIGS. 9(a)-(b) illustrate charts showing measured distributions of spine length and 
volume as a function of time for the population of 52 spines followed in a time-series of 
30 images, and measured distribution of spine motility index fitted to an exponential decay 
function; and 
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FIG. 10 illustrates charts showing lengths of five spines plotted as a function of time 
showing comparison between manual (close circles) and automated (open circles) 
measurements. 

FIG. 1 1(A) depicts three deblurred 2PLSM images of neurons labeled with GFP. Left 
5 image shows neurons cotransfected with wild-type mTOR kinase; center image shows GFP 
transfected control neurons; right image shows neurons cotransfected with an inactive mTOR 
kinase mutant. 

FIG. 1 1(B ) graphically illustrates the results of analysis using an algorithm according 
to the present invention which allows comparison of spine density of GFP transfected neurons, 
10 mTOR transfected neurons and inactive mTOR transfected neurons. 

FIG. 12(A) depicts three deblurred 2PLSM images of neurons labeled with GFP. Left 
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transfected control neurons; right image shows neurons cotransfected with an NLG mutant 
designated AChE. 

15 FIG. 12(B ) graphically illustrates the results of analysis using an algorithm according 

to the present invention which allows comparison of spine density and spine length of GFP 
transfected neurons, NLG transfected neurons and AChE transfected neurons. 

DESCRIPTION OF PREFERRED EMBODIMENTS 

20 The present invention provides efficient, detailed automated analysis of axonic, 

dendritic and spine morphologies, thus allowing assessment of the effects of various genes and 
other substances on neurons. Spines participate in cell-to-cell contact and their growth and 
retraction is powered by actin based motility. The signaling networks underlying spine 
responses are controlled by neurotransmitter receptors and neural activity. Naturally and 

25 artificially induced perturbations in the signaling networks result in alterations in axonic, 
dendritic and spine morphologies and dynamics which can be measured according to the 
present invention with speed and a degree of accuracy and consistency previously unknown. 

An automatic dendritic spine detection and analysis algorithm appropriate for 3-D 
images obtained via laser scanning microscopy or other type of microscopy is presented. The 

30 algorithm of the present invention uses a geometric approach; it is highly automatic and 

contains only a few parameter settings. It can be applied to static images as well as time-series 
images. There is no limitation on the number or the structure of the axons or dendrites in the 
image. In addition to spine length and density, volumetric measurements and spine 



classifications are obtained using the algorithm. Finally, a simple extension of the algorithm 
allows the measurement of distribution of proteins in dendrites and spines. 

The automatic dendritic spine detection and analysis algorithm of the present invention 
offers an objective and consistent analysis, requiring minimal amount of supervision and makes 
5 accessible 3-D morphological characterizations of spine length, volume, shape classification 
and spine density. Comparison of results on spine length and density between the manual and 
automatic approach of the present invention on a large number of samples have validated the 
automatic approach for both static and time-lapse, i.e., time-series, images. 

The automated analysis greatly enhances speed, consistency and objectivity. The timing 
10 results provided below show that automated analysis of time-lapse data consisting of 50 images 
tracking a total of 30-50 spines takes about 4 hours CPU time on a Pentium™ III 500Hz 
processor, whereas for -manual analysis, an experienced user will -typically take 12 to 16 hours. 

For static images the time savings compared to manual analysis may not be as 
significant if experimental conditions vary so significantly that new parameters have to be 
15 determined for each image to be analyzed. However, automated analysis using the algorithm of 
the present invention still provides more detailed, complete and objective quantification than 
manual analysis. 

The automatic algorithm of the present invention assumes that the spines are simply 
connected to the dendrites. Small looping structures in the medial axis indicate pairs of spines 

20 that are too close to be resolved by imaging or segmentation. If it satisfies the protrusion 

criterion each such structures will be detected as a single entity. Resolution of such structures 
as paired spines is currently not implemented. These occurrences are estimated to affect less 
than 2% of the spine population. 

The inventive algorithm enables calculation of spine volume, previously not possible 

25 with manual analysis. Volumetric measurements offer insight into understanding the electrical 
capacity of the spines and the structural and electrophysiological properties of neuronal 
dendrites. It is therefore an important parameter for characterizing dendritic spines. Spine 
volume has not been reported previously in any of the prior art automatic methods. The volume 
measurements agree reasonably well with an SSEM analysis of similarly aged animals, even 

30 though the preparation techniques for the specimens are entirely different. 

The various morphology-based measurement capabilities presented below allows 
application of the investigation of the functional significance of dendritic spines and their 
plasticity to a wide spectrum of experimental and pathological conditions. Automatic 



morphometry significantly improves the scale and accuracy of such studies. Although not 
absolutely essential to the present invention, in certain embodiments, it may be assumed that 
the image to be analyzed is of a bi-phase medium, with one phase being the neuronal cytoplasm 
(dendritic phase), the other being the background tissue. 
5 In accordance with the present invention, a substance is contacted with and/or inserted 

into a neuron and the effects of the substance on morphology is evaluated based on the 
observations relating to spine length and density, volumetric measurements and spine 
classifications according to the present algorithm. After a neuron is analyzed by the present 
algorithm, it may be compared to a control neuron to determine the differences, if any, that are 

10 associated with administration of a substance of interest to the neuron. A control neuron 
includes any neuron used for comparison purposes. For example, a neuron which has been 
subjected to a substance of interest and an imaging dye can be compared to (i) a neuron 
containing that imaging dye or another imaging dye only, (ii) a neuron containing that imaging 
dye or another imaging dye, and another substance of interest, or (iii) a neuron containing that 

1 5 imaging dye or another imaging dye, and a different concentration of the substance of interest. 

Substances which can be investigated include any known material which is amenable to 
be applied to a neuron. Such substances include, but are not limited to, chemical and/or 
physical agents such as microorganisms, DNA, RNA, proteins, peptides, carbohydrates, lipids, 
drugs, radiation, temperature, pH, and diagnostic agents. For example, genes, growth factors, 

20 enzymes, hormones, metals, viruses, bacteria, toxins, dyes, electromagnetism, gamma 

radiation, neurotransmitters, electrolytes, vitamins, minerals, antibiotics, anesthetics, antivirals, 
antiseptics, antimicrobials, antiinflammatories, analgesics, steroids, calcium channel blockers, 
antiarrhythmics, psychotropics, antidepressants and the like may all be contacted with or 
inserted into neurons for evaluation. The substance may be known, or not known, to have 

25 physiological effects since one aim of the present invention is to determine whether or not a 
substance has any effect on a neuron. 

Methods for subjecting neurons to a substance (including causing a substance to 
enter into neuronal cells) are well known. For example, passive or active diffusion utilizing 
concentration gradients (osmolarity), solubilizers, permeation enhancers such as aprotic 

30 solvents, e.g., DMSO, may be utilized to effect entry. 

In a preferred embodiment, genetic engineering techniques are utilized to impart 
substances into neurons. "Transformation" or "transfection" (used interchangeably herein) of 
neuronal cells refers to delivery of nucleic acid (DNA or RNA)) into a neuron by any method. 



Suitable expression vectors include, but are not limited to plasmids, cosmids, phage, 
phagemids, artificial chromosomes and the like. Transfection of host cells can be 
accomplished by, e.g., electroporation, viral transfer, lipid mediated transfer, calcium 
phosphate precipitation, direct injection and biollistic transfer to name a few. 
5 For example, electroporation is suitable for introducing macromolecules, including, but 

not limited to, DNA, RNA, dyes, proteins and other various chemical agents, into neuronal 
cells. Electroporation refers to the permeabilization of cell membranes by application of short 
duration electric field pulses, traditionally between relatively large plate electrodes. During the 
electric pulse, charged macromolecules, including DNA, are actively transported by 

10 electrophoresis across the cell membrane through these pores. Noncharged molecules can also 
enter the pores by passive diffusion. Upon pulse termination, pores reseal over hundreds of 
milliseconds as measured by -recovery of normal membrane conductance values. 

Biollistic transfer is a preferred method for transfection herein and refers to any method 
for introducing foreign molecules into a cell using velocity driven microprojectiles such as 

1 5 tungsten or gold particles. Such velocity driven methods typically originate from pressure 
bursts which include, but are not limited to, helium driven, air driven, and gunpowder driven 
techniques. In biollistic gene transfer, a desired substance, e.g., a plasmid containing a nucleic 
acid sequence of interest, is precipitated on polymeric or metallic beads. Indeed, combinations 
of genes can be delivered by precipitating multiple plasmids onto beads. Thus, for example, 

20 neurons can be transfected with DNA encoding red fluorescent protein to determine cell 

structure, a chimera between a gene of interest and DNA encoding green fluorescent protein to 
track distribution of the gene product and a third plasmid of a specified function. Biollistic 
gene transfer also allows accurate control over the number of plasmids which may be 
introduced into the target cell which allows a measure of control over the amount of target gene 

25 products in the cell. 

Thus, oligonucleotides, chimeric genes, fusion proteins, ligands, receptors, molecular 
labeling systems such as fluorescent molecules, radiolabels, antibodies, antigens, avidin, 
streptavidin, biocytin, and biotin are examples of substances which are suitable for transfection 
herein. Examples of fluorescent molecules (fluorochromes), include green fluorescence protein 

30 (GFP), color shifted mutants of GFP including red shifted mutants, yellow shifted mutants and 
blue shifted mutants, amino coumarin acetic acid (AMCA), fluroscein isothiocyanate (FITC), 
tetramethylchodamine isothiocyanate (TRITC), Texas Red, Cy3.0, Cy5.0 and dextran 
conjugates of fluorochromes. Such labels may be used independently or coupled to other 



molecules such as antibodies, antigens, avidin, streptavidin, and nucleic probes. It is 
contemplated that DNA or RNA encoding fluorescent proteins or other labeling moieties may 
be fused to DNA or RNA encoding a desired protein prior to transfection of such a chimeric 
molecules into the neuronal cell. 
5 The ability to transfer multiple genes into cells is helpful in studies of the interaction 

between different proteins. Transferring genes for colored proteins, such as green fluorescent 
protein, has proven immensely useful for labeling cells in order to visualize their shapes. These 
dyes can also be attached to or otherwise coadministered with other proteins in order to see 
where these proteins are located and move within the cell. As demonstrated below, transfection 
10 was used to fill brain cells with protein dyes in order to observe cell morphology. By 

introducing other genes along with the gene for a fluorescent protein, the effects of various 
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In a preferred embodiment, the neurons being evaluated are studied in preparations that 
are as intact as possible. Accordingly, while individual isolated neurons can be evaluated using 

1 5 the methods described herein, a living brain slice offers an attractive compromise between the 
limitations of cultured neurons (e.g., limited synaptic plasticity) and the experimental 
difficulties of working with intact animals. For example, brain slices may be cultured on 
membranes according to procedures described in Maletic-Savatic et al., Science, 283, pp. 1923- 
1927 (1999) and Stoppinin et al., J. Neurosci. Methods, 37, pp. 173-182 (1991), the contents of 

20 each being incorporated herein by reference. Neurons in cultured brain slices preserve many of 
the aspects of neural functions including normal membrane properties, spine morphologies and 
robust synaptic plasticity. It should be understood, however, that the methods of the present 
invention are also applicable to use in primary cell cultures. 

Any suitable technique involving microscopy to adequately image neurons known to 

25 those skilled in the art may be utilized in accordance with the present invention. Serial section 
election microscopy, confocal laser scanning microscopy (CLSM) and two-photon excitation 
laser scanning microscopy (2PLSM) may be used. A preferred embodiment involves 2PLSM 
since it allows imaging neurons in intact neural tissues. 2PLSM is well-suited to image 
fluorescent molecules. Sectioning may be achieved by limiting the exitation of fluorescence to 

30 a subfemtoliter volume. 2PLSM allows imaging of small fluorescent molecules, e.g., GFP, 

filled structures such as dendritic spines and axon terminals, even in intact brain nervous tissue. 
2PLSM also allows detection of submicron spine length and density changes involved in 
synaptic plasticity in brain slices and the intact brain. Furthermore, 2PLSM allows 

10 



measurement of movement of molecules with submicron resolution. 2PLSM also allows the 
excitation of two fluorophores with different emission wavelengths with the same excitation 
wavelength, facilitating the simultaneous measurement of neuronal structure (in one 
fluorescence channel) and the distribution of selected proteins (detected through a second 
5 fluorescence channel). 

The ability to assess the effects of substances on neurons according to the present 
invention provides an efficient modality to determine what gene products or combinations of 
gene products are involved in the structure of axons and/or dendrites. Indeed, 
structure/function relationships may be uncovered and evaluated according to the present 
10 invention. For example, mutants of naturally occurring proteins can be evaluated for their 

effects on neurons, e.g., dominant negative and constitutively active forms of enzymes. Genes 
that control particular functions in neurons can be identified and studied by transfecting with 
libraries of transgenes and breaking the libraries into smaller pools. Moreover, the effects of 

jQ pharmacologic agents on neurons can evaluated according to the present invention. 
15 The description below describes the algorithm used to analyze 3-D scanning 

^ microscopy images of fluorescent neuronal structures (Section A). The images collected are 

fy analyzed to include, but are not limited to, detection and measurement of dendrites; dendritic 

spines; synapses; distributions of subcellular components using chimeric proteins; and synaptic 

y function using genetically encoded functional probes. The algorithm includes a processing 

M 20 module and an analyzing module for performing the following steps. 

■p. An image is first processed by deconvolution and the dendritic phase extracted (Section 

A.l). The dendrites are identified via their backbones (Section A.2), which are extracted from a 
medial axis construction. Spines are tiny appendages attached to dendrites, and because of their 
small size, they can often only contain a few dye molecules and thus show only dim 

25 fluorescence. They are adjacent to larger, brighter dendrites and they therefore are to detected 
against a hazy background. Unlike the prior art, the spines are not detected from the medial 
axis branches emerging from the backbone as it is difficult to distinguish true spines from 
artifacts by this procedure. Instead, according to the present invention, the spines are detected 
as geometric protrusions relative to the backbone (Section A.3). 

30 Each protrusion is subjected to a protrusion criterion to distinguish true protrusion from 

non-spine dendrite surface irregularity. As very thin necks are too weak to be detectable, some 
spine "heads" appear to be separated from the dendrite and are detected as detached 
components. After initial detection, a search is implemented to associate detached components 

11 



with their appropriate bases. For time-series data, in which the same dendritic branch is imaged 
over a sequence of time intervals, translational effects in time are corrected for and individual 
spines are then traced (Section A.4) through the time ordered sequence of images. Finally, 
morphological characterizations of the population of detected spines are extracted (Section 
A.5). 

The steps described in Sections A.1-A.4 are preferably performed by the processing 
module, while the steps described in Section A.5 are preferably performed by the analyzing 
module of the automatic dendritic spine detection and analysis algorithm 

An imaging, setup and the biological preparation used to obtain data for testing and 
verification of these algorithms are summarized in Section B and the Examples herein. Results 
of the application of these algorithms to the analysis of hippocampal CA1 neurons and a small 
number ofhippocampalC A3 neurons are presented in the Examples as well. 
A. Image Analysis 

A.l Image Decon volution and Segmentation 

The intrinsic spatial resolution limits of optical microscopy arise from the diffraction of 
light; light from a point source is ideally imaged to a larger spot characterized by the Airy 
function. The measured spread resulting from a given optical setup is referred to as the point- 
spread-function (PSF). As a result, the intensity recorded in any voxel (volume element) of a 
digitized image is a convolution of intensities from its neighborhood. 

Deconvolution is used to correct aspects of the image degradation due to the PSF. A 
variety of deconvolution techniques are available which employ either theoretical or 
experimental measures of the PSF. In addition, blind deconvolution methods can be employed 
which, concurrently with the deconvolution, reconstruct an estimated PSF of the image. 
However, the presence of noise and the band-limited nature of the PSF limit the improvements 
by means of classical deconvolution techniques. Therefore some blurring will remain even after 
deconvolution due to a trade-off between sharpening of the image and noise amplification. 

In addition, the photomultiplier tube (PMT) detectors used in most laser scanning 
microscopes are "noisy". Even in darkness PMTs produce spontaneous bright pixels which 
may be referred to as "shot noise". One can deal with shot noise by applying a median filter to 
the image. Median filters are known in the art. See, e.g., Tukey, Exploratory Data Analysis, 
Addison- Wesley, Reading, MA (1971). For example, the median filter may be a non-linear, 
lowpass filter which replaces the greyscale value of each voxel v in the digitized image by the 
median greyscale value of v and its 26 neighbors. This effectively removes shot noise but not 

12 



real spines which, under typical magnifications employed, have an effective width covering 
many voxels. 

A prior art iterative reblurring deconvolution algorithm was applied to the median 
filtered image, which requires either a theoretical or experimentally measured PSF. See, 
5 Kawata and Ichioka, J. Opt. Soc. Am., 1980; 70:762-772, the contents of which are 
incorporated herein by reference. Briefly, iterative reblurring proceeds as follows. Let 

o (0) (x,y y z) denote the experimental 3-D image and h (x y y,z) an appropriate PSF. Let * and n 

denote the convolution and correlation operators. The deconvolved image 5 {k) (x,y,z) in the k- 
th iteration is 

10 5 {k) = <5 ( *-° + { o (0) n h - <3 ( *-° *(hnh)}. 

A non-negativity constraint is applied to 6 (k) at the end of each iteration. 

For the images analyzed in the Examples herein, the PSF was measured by imaging a 
Q number of sub-resolution microspheres and averaging their individual PSFs to reduce noise. 

FIG. 1 demonstrates the result of iterative reblurring a raw image (FIG. 1(a)) after 5 (FIG. 1(b)) 
-i; 15 and 20 (FIG. 1(c)) iterations. Preferably, one employs k^ = 5 iterations of deblurring. It 

should be understood that those skilled in the art may use other deconvolution algorithms for 
gQ reblurring. 

Segmentation is a generic imaging term for labeling each voxel in a greyscale (or color) 
*F image with an integer identifier designating its "population type". For dendritic morphometry, 
p 20 this requires distinguishing neuron voxels from the background tissue voxels. A large number 
f: of segmentation algorithms are available. See, e.g., Pal and Pal, Pattern Recognition, 1993; 
26; 1277-1294. As the dendritic images are processed first by median filtering and 
deconvolution, simple thresholding is used for this final segmentation step; all voxels of 
intensity greater than a threshold value are identified as neuron, otherwise as background. In 
25 general, a trade-off between selection of dim spines and reduction of noise on the dendrite 
surface is made in selecting the threshold. 
A.2 Dendritic Backbone Extraction 

Geometric analysis of a 3-D irregularly shaped object is difficult; such analyses 
typically employ models based upon geometrically simple "unit" objects. The algorithm of the 
30 present invention builds upon the medial axis algorithm disclosed by Lee et aL 9 "Building 
skeleton models via 3-D medial surface/axis thinning algorithms," CVGIP: Graph. Models 
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Image Process., (1994); vol. 56: 462-478, to provide a skeleton from which the backbone of 
each dendrite in an image can be extracted. 

Intuitively, the medial axis captures a geometrically faithful skeleton (consisting of 
curve segments joining at vertices) of an object. In a digitized image these curve segments 
5 consist of linked sequences of voxels, with the vertices being voxels at which these segments 
join together. An example of the medial axis of a portion of a segmented dendritic image is 
shown in FIG. 2(a). (This is a view perpendicular to the optical axis.) The medial axis obtained 
for the dendritic phase contains the backbone ("centerline") of each dendrite as a subset. 
In addition to the backbone however, the medial axis contains "spurs" and other 
10 features which correspond to spine-related or non-spine-related surface features (e.g., incipient 
dendritic branches); to surface artifacts resulting from digitization effects, segmentation errors, 
and boundary effects due to the finite imaged volume; or to spurious cell "debris". Due to 
^ resolution limits, spines emerging near each other may appear to have overlapping tips in the 
y3 digitized image, resulting in the appearance of small loops in the medial axis (top of the parent 
15 branch in FIG. 2(a)). A separate skeleton for each disconnected component of dendritic phase is 
■0 also contained in the medial axis. 

fy From the medial axis, the backbone for each dendrite is extracted (FIG. 2(b)) in two 

,y=ps steps. The first step is achieved by removing the medial axis segments corresponding to all 
.0 disconnected dendritic components, and trimming short spurs and loops on the medial axis. 
^ 20 Removal of long "spurs" is problematic as they may correspond to filopodia or incipient 
y branches of the dendrite. These long "spurs" are dealt with in the next step. 
;M In the second step, a backbone is traced through each dendritic branch employing a 

decision based upon minimum deviation angle whenever a vertex on the trimmed skeleton is 
encountered. If necessary, the number, n, of dendrites in the image can be specified so that only 
25 the n longest backbones are retained. Any remaining medial axis segment that is not part of a 
traced backbone removed. The final set of dendritic backbones extracted from FIG. 2(a) is 
shown in FIG. 2(b). 
A.3 Spine Detection 

With the dendritic backbones isolated, the spine detection algorithm of the present 
30 invention proceeds in four steps: detection of detached dendritic phase components (Section 
A.3. a); detection of attached spine components (Section A.3.b); elimination of spurious or 
incomplete spine components (Section A..3.c); and merging of spine components (Section 
A.3.d). Since a spine may be composed of one or more detached pieces and possibly an 
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attached base in the segmented image, the identification of any spine is not finalized until all 

four steps have been completed. 

A.3.a Detached spine component detection 

Dendritic phase components disconnected from the backbone-containing dendrites are 
detected and tentatively identified as detached spine components. For each detached spine 
component a record is kept of its center of mass, the closest dendritic backbone voxel and the 
dendritic surface voxel lying on the line joining the center of mass to the backbone voxel. 
Detached dendritic phase components that are further from the nearest dendrite surface voxel 
than a maximum distance tolerance are interpreted as false positive signals and ignored. For 
the images analyzed in the Examples herein, the length tolerance is 6 pm. It should be 
understood that the length tolerance can be adjusted when appropriate. 
A.3.b Attached spine component detection 

Ignoring the detached spine components, every dendrite phase voxel v is labeled with a 
distance, <4(v), to its closest backbone voxel. Thus tips of protrusions on the dendrite surface 
are assigned the largest distances. These tip voxel locations are then processed in descending 
order of 

For each tip voxel S , a sequence {C,} 5 i = 1,..., d b (s), of candidate spines is generated. 
Candidate C i consists of all voxels w whose distance d s (w) from S is d s (w) < i . Fig. 3(a) 
shows a 2-D projected view of the candidates (gray voxels) C 13 ..., C n for a tip S having d b = 
12. The smaller candidates clearly contain insufficient voxels to correctly represent the spine, 
whereas the larger candidates protrude too far below the dendrite surface. The optimal choice 
of a spine candidate would terminate at the surface of the dendrite. This is achieved by 
estimating for the local thickness of the dendrite as explained below. 

To estimate the local thickness and choose the optimal candidate, a ring of "spine- 
surface boundary points" are determined for each candidate. For clarity of explanation, the 
algorithm is first illustrated in 2-D assuming the image is projected onto its focal plane. The 3- 
D algorithm is described afterwards. Two "surface boundary voxels" P x and P 2 on the 2-D 
projection are shown in FIG. 3(b). For each candidate C, , along with the two surface boundary 
voxels P x c ' and P 2 q , a "base voxel" E q having the furthest penetration into the dendrite is 
also determined. An idealized sketch of such a projection illustrating S 9 P x , P 2 and E is 
shown in FIG. 3(b). These surface points are used to determine the best measure of the dendrite 
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thickness as follows. A "reference candidate" C R is selected to be the smallest volume 

candidate with the minimum surface to backbone distance, d p = min,. { d b {pf 1 ), d b (p 2 C/ ) } . 

A reference candidate C R and the distance d p are illustrated in FIG. 3(b). 

In an ideal case 2 d p is the width of the dendrite and the best candidate for the spine 

5 would be C i where j = d b (s) - d p + 1 . In practice, spines and dendrite surfaces in the images 

are quite irregular and additional criteria must be satisfied before a final candidate is accepted. 
To be accepted as a true spine, the candidate is required to satisfy a heuristic protrusion 

criterion. LetD^ pp denote the perpendicular distance from S to the line segment Pf'Pf 1 of 

C, and let D^ pp denote the perpendicular distance from the base voxel E Ct to the line 

10 segment Pf'P^ as illustrated in FIG. 3(c). The spine candidate C. is required to satisfy 

DsUw * DeUp A (1) 
Criterion (1) clearly requires that favorable spine candidates protrude out further from 
the dendrite surface than they protrude into it. It however favors spines with narrow necks, as 
the following argument shows. Consider an ideal, symmetric spine whose base is the arc of a 
15 circle with radius R as sketched in FIG. 3(c). (Clearly i plays the integer role of the radius for 
the digitized candidates C i .) 

Let 2 W denote the distance between P x and P 2 . Simple trigonometry gives 

D E ^ PA =R -JR 2 -W 2 , (2) 

D s _ PA = ylR 2 -W 2 . (3) 

20 Then the protrusion criterion (1) requires R > 2W I in order to accept a spine candidate, 
favoring narrow-necked spines and rejecting wide-necked spines. Since R will increase faster 

than Wfor true spines (i.e., i increases faster than the distance P x Ci P 2 C/ ), instead of testing only 
the candidate C } for the protrusion criterion, a range of candidates C i are tested for which i is 
close to y. Specifically, candidates having values i in the range 
25 d b (S) -d p <i<d b (S) - (d p + d e - \)l 2 , where d e =d b (E c "), are considered. If one of the 

candidates C i with i value in this range satisfies (1), then the candidate C . is accepted as a true 
spine; otherwise the protrusion is rejected as a spine. 
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The 3-D algorithm proceeds in a similar fashion except that instead of using the 
projected pair of surface boundary points, the entire ring of surface boundary points in 3-D are 
considered. The minimum "surface points to backbone" distance is used to find the reference 
candidate C R to correct for the local thickness of the dendrite. The distance from S (or E) to 
5 the ring of voxels is calculated by measuring the perpendicular distance of S (or E ) to the 
plane that best fits the ring of voxels in 3-D. 
A.3.C Component Elimination 

Spines touching the boundary of the imaged region are ignored as they are incomplete. 
This is also a useful technique for eliminating "debris" and other axons or dendrites in the 

10 background of the image that are near or touching the dendrites of interest. The inventive 

algorithms are written to allow for the imposition one or more non-overlapping polygonal areas 
on the plane of the image slices. The interior of the union of these polygons is regarded as the 
region of interest for the spine detection algorithm; any structure exterior to the polygons is 
ignored. By setting the polygonal edges to cross through unwanted structures they are also 

15 automatically ignored. As mentioned, detached components further from the dendrite surface 
than a maximum distance are also eliminated. 
A.3.d Component Merging 

As a spine may be identified from multiple detached "head" and attached "base" 
components, a final merging algorithm which accounts for the position and orientation of all 

20 possible spine pieces is performed. The merging algorithm considers every component, 

checking for possible merges with other components. Any merged entity is reconsidered as a 
new single component, and rechecked for possible further merges. 

Merging can occur between two detached (DD) components (the merged entity is still 
considered a detached component) or between detached and attached (DA) components (the 

25 merged entity is then considered to be an attached component). Two criteria are employed for 
DD or DA type merging. 

The first criterion is maximum separation; the two components to be merged are 
required to be close enough (a center-of-mass to center-of-mass separation < 3|im). The second 
criterion requires appropriate relative orientation of the two components as demonstrated in 2- 

30 D in FIG. 4. For DA type merging, the tip S of the attached component A is required to lie 

within the triangle DP\P 2 . (In 3-D, the tip S is required to lie within the cone determined by D 
and the ring of spine-surface boundary points.) For DD type merging, the average angle 



17 



5? 15 



subtended by the center of mass of each spine with the surface voxel locations of both spines is 
required to be less than 30°. 
A.4 Image Registration and Spine Tracing 

A time sequence of 3-D images must be registered to correct for possible translational 
5 movement of the specimen. After registration, individual spines are then traced and identified 
through the image sequence. 

Each consecutive pair of images F i and F i+l are co-registered using the spines 
separately identified in each image. The offset 5 = [o x ,o y9 o 2 ) of F M with respect to F i is 

allowed to vary within a window \o x \ < w x , jo^ | < w y , \o 2 \ < w 2 . Only integer voxel offsets are 

10 considered. A conventional registration method maximizes the cross correlation of two images; 
thus no decision can be made until the correlation arrays are computed for all offsets. Instead, 
an efficient sequential search method, as disclosed by Barnea and Silverman, 1972 in "A class 
of algorithms for fast image registration," IEEE Trans. Computers, 1972, vol. C-21, 179-186, is 
utilized which computes the l x norm (absolute value sum) image difference 

fi over all offsets o in the window for which e(5 ) is less than a predetermined threshold value T. 

2 The offset 5 with minimum e(o ) provides the optimal registration. In practice, 

*7 w x = w y = w 2 = 5 voxels, and Tis the average number of total spine voxels in F { and F i+l . 

if Individual spines are traced through the time-series. Two spines at different times are 

20 considered to be the same if their percentage overlap (measured in voxels) is larger than 25% of 
the volume of at least one of them. 
A.5 Morphological Characterization 

Spine length, density and volume are computed in accordance with one aspect of the 
present invention. Spines are also classified according to their shape. 
25 For a detached spine (without any attached component), the spine length is determined 

by the distance from the recorded dendrite surface voxel (corresponding to the associated 
dendrite) to the furthest spine voxel (corresponding to the detached spine) from the dendrite. 
For spines that are fully or partially attached (consisting of a base and one or more detached 
components) to the dendrite, the spine length is determined by the distance from the center of 
30 mass of the base boundary points to the furthest spine voxel (possibly detached from the 
dendrite). For the images analyzed in the Examples herein, the automated spine length 
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measurement is calculated from a 2-D projection. The reason for this is that the manual spine 
analysis measurements against which the automatic analysis results are to be compared are 
performed in 2-D by projecting the 3-D stack of image slices along the optical direction. 
Spine density is computed as the number of spines per unit length of dendritic 
5 backbone. For purposes of comparison with manually analyzed images which are analyzed in 
2-D projection only, backbone length is also measured from a 2-D projection onto the slice 
plane. 

Spine volume is measured according to the intensity values of the deconvolved 
greyscale image. For 2PLSM, the excitation of fluorescence is limited to a sub-femtoliter focal 

10 volume ( ~ 0.5 x 0.5 x 1 .5 jam 3 ) which is larger than that of individual spines. The intensity 
value recorded for each voxel in a spine is a sum of the fluorescence from all dye molecules 
excited within the focal volume. The maximum intensity voxel near the center of a spine is 
therefore a measure of the volume of a spine. As the larger cross-sectional areas of a dendrite 
are typically larger than the maximum cross-sectional area of the focal region, the maximum 

1 5 voxel intensity recorded along the dendrite backbone is a measure of the size of the focal 

volume, assuming the fluorescence is saturated near the center of the dendrite. See, Svoboda et 
al., Science, 272, pp. 589-593 (1996) and Sabatini and Svoboda, Nature, 408, pp. 589-593 
(2000). The spine volume is defined as the ratio of the maximum spine intensity to the 
maximum dendrite intensity multiplied by an empirically determined focal volume, 

™ o . t, t Maximum Spine Intensity _ f __ _ 

20 Spine Volume = x Focal Volume . 

Maximum Dendrite Intensity 

The following classification of spine shapes is used stubby, thin, mushroom. Spine shape is 

decided based on spine length (Z), head diameter (d h ) and neck diameter (d n ). In general 

terms for thin spines, spine length should be much greater than the neck diameter (l » d n ) . 
For mushroom spines, spine length should not exceed neck diameter by more than a factor of 5, 
25 and the head diameter should be much greater than the neck diameter (d h »d n ). For stubby 
spines, the neck diameter is approximately equal to the length of the spine. The specific criteria 
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adopted in this classification utilize the ratios Lld n and d h ld n to classify theirshape as 
summarized in Table 1 . 

Table 1: Ratio criteria for the classification of stubby, thin and mushroom spines. 
5 ~Lld n 7jd n 





[0,1.3) 


[13,3) 




[0,2/3) 


stubby 


mushroom 


mushroom 


[2/3,2) 


stubby 


stubby 


stubby 


[2,3) 


stubby 


mushroom 


mushroom 


[3,5) 


thin 


mushroom 


mushroom 


[5,») 


thin 


thin 


thin 
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A.6 Measurement of the distribution of fluorescent proteins. 

q Segmentation is performed on data from one fluorescence channel. The other 

yQ fluorescence channel can be used to measure the distribution of protein components, such as 

ffl 20 chimeric proteins linked to GFP. The analysis is simply to measure the fluorescence on the 

fJ second channel as a measure of protein concentration, in pixels that were previously 

S determined to belong to the neuronal structure. 

m 

' B, Image Acquisition 

a 25 From a data analysis standpoint CLSM and 2PLSM provide essentially equivalent 

]* challenges. However, to gain an understanding of the dynamics of neuronal circuits, neurons 
should preferably be studied in preparations that are as intact as possible. For many questions 
J3 of sub-cellular physiology, as stated above, the living brain slice offers an attractive 
" compromise between the obvious limitations of cultured dissociated neurons and the 
30 experimental difficulties encountered when working with intact animals. 

One problem with brain slice physiology has been that scattering of light makes 
traditional optical microscopies, including CLSM, difficult in living tissues. For these reasons, 
as also stated above, the data is preferably collected using 2PLSM, which allows high 
resolution fluorescence imaging in brain slices up to several hundred microns deep with 
35 minimal photodamage. 

At present, parameters that require routine adjustments include the region of interest 
and segmentation threshold. Other parameters that are used in the deconvolution, backbone 
extraction, spine component elimination and tracing algorithms are empirically determined; 
they remained the same throughout all of the examples that have been described herein. 
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Segmentation is crucial to the analysis due to the relatively low intensity associated with small 
spines and the high intensity of the dendrites. Choosing a critical threshold is important; simple 
thresholding is adequate for most images that have been preprocessed by median filtering and 
deconvolution. 

5 It is contemplated to provide the algorithm of the present invention within a server 

accessible via a network, such as the Internet or a local area network (LAN), by a plurality of 
users. The users can then transmit data to the server for processing using the algorithm. The 
results would then be transmitted by the server to the user via e-mail or by other known 
techniques. 

10 The following examples are included for purposes of illustrating certain aspects of the 

present invention and are not intended to limit the invention as defined by the claims herein. 

EXAMPLE 1 



m is 



Sample Preparation and Microscopy 



Cultured hippocampal brain slices were prepared from 7 day old rats. After five days in 
vitro a small subset of neurons were biollistically transfected (for example, Lo, et al., (1994). 
Neuronal transfection in brain slices using particle-mediated gene transfer, Neuron 13, 1263- 
1268, the contents of which are hereby incorporated by reference) with a plasmid carrying the 
g 20 gene for enhanced green fluorescent protein (GFP) (commercially available from Clonetech) At 
least two days after transfection, slices were transferred to a perfusion chamber for imaging. 
Labeled neurons were identified and imaged using a custom-made 2PLSM laser scanning 
microscope (as described in Mainen, Z. F., Maletic-Savatic, M., Shi, S. H., Hayashi, Y., 
Malinow, R., and Svoboda, K. (1999), Two-photon imaging in living brain slices, Methods 75, 
25 231-239, the contents of which are hereby incorporated by reference). The light source was a 
Ti:sapphire laser running at a wavelength of «990 nm (repetition frequency 80 MHz; pulse 
length 150 fs). The average power delivered to the backfocal plane of the objective (40x, NA 
0.8) varied depending on the imaging depth (range 30 to 150 mW). Fluorescence was detected 
in whole-field detection mode with a photomultiplier tube. 
30 EXAMPLE 2 

Static Analysis 

To validate the automatic spine detection algorithm, an experiment, E\ 9 involving -200 
spine measurements over 15 dendrites of hippocampal CA1 and a small number of CA3 



21 



neurons was performed. The same imaged regions were subjected to both automatic and 
manual analysis. A total of 174 spines were identified by both methods; an additional 10 spines 
were identified only by the manual method; and a further 28 spines were identified only by the 
automatic method. 

5 The results of the manual and automated analysis for one of the dendrites in this 

experiment are illustrated in FIG. 5(a)-(b). Twenty-one spines were detected by both methods; 
three additional, relatively short, spines were detected by the automatic method. FIG. 6 
compares the individual spine lengths, average spine length and spine density measured for this 
particular dendrite. Spines 8 and 22 demonstrate the difficulties encountered by both detection 
10 methods when two spines appear to overlap. For spine 8, the manual detection has identified 
only the shorter of two spines which appear to overlap; whereas the automatic method has 

iH*»«tifi<aH r»nl\/ tVi^» 1rm nf±r T*nr er\ir»p> 99 afrcnn Hn/n cninpc annpar tn nvprlan jmH arp r.nn QiHprpH 

as a single spine by the automatic method. On the other hand, the manual method failed to 
identify either of them. 

15 Table 2: Measured mean spine lengths ( ± standard deviation) for the spines in experiment E\. 

Population size Method Mean spine length 

(urn) 

20 174 Manual 1.05 ± 0.62 

174 Automatic 1.08 ± 0.63 

10 Manual 0.75 ± 0.57 

28 Automatic 0.38 ± 0.28 

25 Table 2 compares the mean spine length measured by each method for the population of 

174 spines detected in common. The mean lengths for those spine detected by only one of the 
methods are also presented. For the common detected spine population, the manual and 
automatic spine length measurements agree to within one standard deviation, though the 
standard deviations are large. (The large standard deviation is partly due to averaging over 

30 spines of different shape classification.) A paired samples Student's t-test to determine whether 
the difference in measurements by the two methods is significant provides a stronger test of the 
agreement between the two methods of length measurement. Column 1 of Table 3 summarizes 
the results of the paired t-test; there is no significant difference between the two methods of 
length measurement. 

35 
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Table 3: Paired samples t-test results for measured spine lengths and densities, experiment E\. 

Spine length Spine density 



5 Degrees of freedom 173 14 

t statistic -1.24 -0.87 

p-value 0.22, two-sided 0.40, two-sided 

A one-way ANOVA was used to test for any dependence of dendrite origin on the 
10 observed differences in measured spine length. The test produces an F statistic value of 0.95 

(d num = 14 and d den = 159) with a p-value of 0.51. Thus the differences in spine length 
measurements between the two methods are uniform across the different dendrites. 

The mean spine length for the 28 spines detected only by the automatic method of the 
present invention indicates a population of smaller length spines. The results in Table 4 of an 
15 independent samples t-test for these 28 spines shows that the difference in spine lengths is 

significant compared to that obtained by either measurement method for the population of 174 
commonly detected spines. The results indicate that the automatic algorithm of the present 
invention is detecting short spines more consistently than the manual method. 
Table 4: Independent samples t-test results, experiment E\. 

20 



Method (population 


Manual (n = 174) 


Automatic (n = 174) 


Automatic only (n = 28) 


df = 200 


df= 


= 2 




t = 5.63 


t = 


5.82 


P = 


6 x 10" 8 , two-sided 


P = 


2 x 10" 8 , two-sided 


Manual only (« = 10) 


df = 182 


df= 


= 182 




t = 1.51 


t = 


1.65 


P = 


0.13, two-sided 


P = 


0.10, two-sided 



For the 10 spines detected only by the manual method, the mean length measurement 
lies midway between that obtained for the common and automatic only populations. An 
independent samples t-test (Table 4) shows no significant differences with the measurements 
35 obtained by either method for the 174 commonly detected spines. The results (df = 36, t = 2.65, 
p = 0.01, 2-sided) of an independent samples t-test between these 10 spines and the 28 detected 
only by the automatic algorithm indicate a significant difference between the spine lengths of 
these two populations. 
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Visuaf observation of these 10 spines reveals that 7 of them touched the boundary of the 
image region and were consequently rejected by the automatic algorithm. The remaining 3 
were not resolved by the automatic algorithm as each touched some neighboring spine (which 
was detected). This is thus a reflection of the combined effectiveness of the median filter, 
5 deconvolution and simple thresholding algorithms in segmenting the images. Based upon visual 
investigation of the images, no more than 2% of the spines were estimated not to have been 
resolved by the automatic algorithm due to segmentation related effects. 

Table 5: Measured mean spine density (± standard deviation for the dendrites in experiment 
10 E u 

Population size Method Mean spine density 



15 15 Manual 0.45 ± 0.09 

15 Automatic 0.47 ±0.15 

Table 5 compares the mean spine densities separately measured by each method for the 
15 dendrites. (For the manual method this is a total population of 184 spines; for the automatic 

20 method, 202 spines.) The density measurements by either method agree to within one standard 
deviation. For the paired sample of 15 dendrites, a Kolmogorov-Smirnov test shows that the 
dendrite by dendrite difference in the automatic and manual measured densities is very close to 
normal, so that a paired-dendrite samples t-test can be applied. Column 2 of Table 3 
summarizes the result; the automatic spine density measurement is not significantly different 

25 from the manual. 

EXAMPLE 3 
Static Analysis 

30 For spine volume measurement and shape classification, automated results are reported 

herein as no manual determination is available. A second experiment, Ei, was performed under 
the same experimental conditions to increase the sample size (E\ + Ei) up to -700 spines. The 
spine volumes were calculated from the ratio of the maximum intensity values of the spine to 
the dendrite, as described in Section A.5 using an empirically determined focal volume of 0.5 x 

35 0.5 x 1.5 (xm 3 . FIG. 7 shows the volume-length correlation plot for the spines measured in 
experiments E\ and E2 according to their determined classification. The mushroom shaped 
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spines occupy the widest spectrum of lengths and volumes. The ratio of stubby : mushroom : 



thin spines is 0.54 : 0.36 : 0.10. 

Table 6: Comparison of automated spine length and volume measurements of hippocampus 
cells in PND 7 cultured neurons with the PND 15 SSEM results of Harris et al y "Three- 
5 dimensional structure of dendritic spines and synapses in rat hippocampus (CA1) as postnatal 
day 15 and adult ages: implications for the maturation of synaptic physiology and long-term 
potentiation," J. Neurosci. 1992, vol. 12, 2685-2705. 



Measurement 


Method 


Shape classification 








stubby 


mushroom 


thin 


volume (um 3 ) 


automatic 


0.07 + 0.04 


0.06 ± 0.05 


0.06 ± 0.04 




SSEM 


0.11 ± 0.07 


0.18 ± 0.09 


0.05 ± 0.03 


length (um) 


automatic 


0.65 ± 0.37 


1.35 ± 0.55 


1.38 ± 0.54 




SSEM 


0.65 + 0.38 


0.95 ± 0.30 


1.40 ± 0.39 



Table 6 summarizes the average spine volume and length measurements in each shape 
20 category and presents comparison with the SSEM results (Harris et ai 1992; Table 4) on rat 
hippocampus CA1 cells for postnatal day (PND) 15 animals. All automatic measurements are 
within 1.5 standard deviations of the SSEM result, though the volume results are generally 
smaller. It is noted however that the automatic results come from cultured neurons and younger 
aged animals. In addition, no corrections for any fixation-induced changes have been 
25 performed in the SSEM study. 

EXAMPLE 4 
Time series analysis 

Time-series data provides the ability to capture dynamic changes in dendritic spine 
30 morphology. A series of 50 3-D images was taken at 30 second intervals spanning a time 
period of 25 minutes. FIG. 8(a) shows the number of spines detected in the images as a 
function of time using the automatic method. On average 27.5 ± 2.3 spines were detected in 
each image. An interest lies in the question of the frequency of observation of any particular 
spine over time. 

35 In total, 52 spines were detected and traced through the time-series. FIG. 8(b) shows 

how the observations of the 52 spines were distributed in time. The spines were indexed 
(1 52) according to the first time in which they appeared. Thus, 27 spines were observed over 
the full 25 minutes of image taking; among those, 16 were present at all time points. 
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FIG. 9(a) summarizes the distribution of spine length and volume as a function of time. 
Both length and volume distributions are skewed, with smaller lengths and volumes 
dominating, consistent with the stubby, mushroom, and thin spine ratios noted above. The 
dynamics of the spines are measured using an index for spine motility, which is defined as the 
5 summed difference in length of a spine in time divided by the total number of time steps. 
FIG. 9(b) plots the distribution of motilities for the 52 spines traced in this series. For this 
limited data set, the number of spines ( n ) decreases with motility (m ) approximately as 

n{m)=n{0)e-> 69m . 

A comparison between automated and manual length measurements was made for a 
10 limited subset of this time series data; manual length measurements were made on a subset of 
10 of the 52 spines. FIG. 10 presents comparisons between the automated and manual spine 
length measurements as a function of time for 5 of the spines chosen to represent different 
f=^ average lengths. Consistently longer lengths were measured by the manual method for the 
!fy longer spines (1 and 2). For the medium length spines (3 and 4), the manual and automated 
fU 15 results are very similar. For the short spine (5), some deviations are observed; occasionally the 
^ spines were not detected by the automatic method. 

fy A paired t-test was performed on this subset of 10 spines to determine whether the 

- average spine length determined by the automated measurement is significantly different from 

lb 2 that determined by the manual measurement. A Kolmogorov-Smirnov test shows that this 
;H= 20 sample of 10 difference measurements (-0.32 ±0.17 |im) is very close to normal so that a 
□ paired samples t-test can be applied. The observed t statistic value is -2.79, with p = 0.02, two- 
;F= sided, revealing some significance in the averaged difference. 

While this is contrary to the results in experiment E\ 9 it is noted that the manual 
measurements were made by a different user than in Ej. Pearson correlation was therefore used 
25 to test whether the automatic and manual measurements are correlated in time. The correlation 
values (r) obtained for the 10 spines for n = 50 time points range from 0.92 to 0.29; with two- 
sided p-values ranging from 0.00 to 0.04. Thus, for these 10 spines significant correlations 
between the automatic measurements and the manual measurements in time are observed. 
Therefore, the existence of a systematic bias between the manual and automated measurements 
30 for this set of data that can be attributed to a change in the user making the manual 
measurements. The significant Pearson correlation however indicates that the manual 
measurements are duplicating the trends found by the automatic measurements. 
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EXAMPLE 5 



Evaluation of the Effects of mTOR Kinase on Spine Formation 

mTOR kinase controls the phosphorylation of the translation regulators p70 56k and 4E- 
5 BP 1 . It is a central regulator of cell growth and highly expressed in dendrites and highly 

associated with synaptic proteins. Rat hippocampal CA1 neurons were transfected by biollistic 
gene transfer with GFP alone or together with the wild-type kinase (mTOR wt) or a non- 
functional mutant (mTOR kd) that acts as a dominant negative (for details of plasmids see: 
Sabatini et al., (1999). Interaction of RAFT1 with gephyrin required for rapamycin-sensitive 
10 signaling, Science 284, 1161-4, the contents of which are hereby incorporated by reference) 
(for methods see Example 1). The 3d structure was imaged using 2PLSM as described in 
Example 1 and the inventive algorithm utilized to analyze dendritric shape. As can be seen in 
FIGS. 1 1(A) and 1 1(B) mTOR kinase controls spine size and density, i.e., mTOR wt increased 
spine density and mTOR kd decreased spine density as compared to the GFP control. Data 



CD 

»1 15 from at least 5 neurons were in every group. 



ry EXAMPLE 6 

^ Evaluation of the Effects of Overexpression 

□ of Neuroligin on Spine Formation 

u 20 Neuroligin (NLG) is believed to be involved in synapse formation. When nonneuronal 

•tf cells are engineered and cultured to express NLG, they are associated with the development of 
H presynaptic structures in contacting axons, suggesting that NLG-neurexin interactions are a key 
step in synapse formation. Rat hippocampal CA1 neurons were transfected by biollistic gene 
transfer with GFP alone or together with NLG or a mutant form of NLG that does not contain 
25 the AchE domain (designated "AChE") (See, Scheiffele et al., (2000). Neuroligin expressed in 
nonneuronal cells triggers presynaptic development in contacting axons, Cell 101, 657-69, the 
contents of which are hereby incorporated by reference). The 3d structure was imaged using 
2PLSM as described in Example 1 and the inventive algorithm utilized to analyze dendritic 
shape. As can be seen in FIGS. 12(A) and 12(B), overexpression of NLG did not produce a 
30 detectable phenotype change in spine density or spine shape. Data from at least 5 neurons are 
in every group. 
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EXAMPLE 7 

Measurement of the distribution of GluRl-GFP in dendritic spines. 

Glurl is a synaptic receptor that is thought to play an important role in synaptic plasticity. To 
measure the distribution of Glurl in dendritic spines neurons were transfected with a virus 
5 expressing Glurl -GFP (for methods see Shi, S. H., Hayashi, Y., Petralia, R. S., Zaman, S. H., 
Wenthold, R. J., Svoboda, K., and Malinow, R. (1999). Rapid Spine Delivery and 
Redistribution of AMPA Receptors After Synaptic NMDA Receptor Activation, Science 284, 
1811-1816, the contents of which are hereby incorporated by reference). Neurons expressing 
glurl -GFP were then patch-clamped and filled with a red fluorohore (Texas Red, Molecular 
10 Probes). A two color image was then acquired using 2PLSM, with an excitation wavelength of 
910 nm. The red image was used for segmentation and to perform a spine analysis. The green 
image was used to estimate the distribution of Glurl -GFP in dendrites and spines. 

It will be understood that various modifications may be made to the embodiments and 

jrsl 

jl examples disclosed herein. For example, alternative algorithms can be created to accomplish 
SI 15 the criteria set forth above. It should be understood that, notwithstanding the emphasis on 
U spine morphology, the algorithms described herein can be applied to the larger axonic, 
pj dendritic and cell body structures of neurons to determine length, volume, shape classification 
^ and density. Therefore, the above description should not be viewed as limiting, but merely as 
□ exemplifications of preferred embodiments. Those skilled in the art will envision other 
2 20 modifications within the scope and spirit of the claims appended hereto. 
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